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Bosonic Structure of a 2-Dimensional Fermion Model 
with Interaction among Different Speices 



Jiro SakamotoO and Yasunari Heike 
Physics department, Shimane University, Matsue 690-8504, Jo-pan 

We study a massive Thirring-like model in 2-dimensional space-time, which contains two 
different species of fermions. This model is a field theoretical version of the quantum me- 
■ chanical model originally proposed by Glockle, Nogami and Fukui, where different fermions 

interact with each other through 5-function potentials. We derive a corresponding boson 
model by the bosonization technique in the path integral formulation. This is a simple but 
non-trivial extension of the freedom of the bosonization technique. Operator correspon- 
dences between fermion and boson fields are given. One of these could not be realistically 
expected from the naive correspondence of the original single-species models. It is essential 
for this point that in our model fermions of the same kind do not interact with each other 
directly. We find that for a specific value of the coupling constant, one boson field becomes 
free while the other is a Sine-Gordon field. For this case, therefore, our two-species model is 
equivalent to the ordinary Sine-Gordon model of a single boson field. 



> 

CO , 

■ §1. Introduction 

In previous papers we have studied relativistic bound states of a 1-space quan- 
tum mechanical system containing different species of massive fermions in order to 
investigate the relativistic effects for such composite systems. EP This model is an ex- 
tended one from the model of two kinds of fermions originally proposed by Glockle, 
O ■ Nogami and Fukui (GNF).& The Hamiltonian of this model is given by 

H = ^{-ictiPi + m(3i} - - ^(1 - a-joy)^ - Xj), (M 

i=l i^j 

where i and j denote fermion species. It is essential in this model that different 
fermions interact with each other through the 5- function potentials, while fermions 
of the same kind do not interact with each other directly. All the requirements of 
quantum mechanics and special relativity are satisfied. We found an exact solution 
for n-body bound state which contains n different particles. 

The GNF model is, however, based on the single-electron theory, where anti- 
particles are not supposed to exist, and necessarily its Hamiltonian (IT) is not 



positive-definite. One way to overcome this defect would be to go into field theory. 
It is seen that the GNF model can be derived from a massive Thirring-like model in 
2-dimensional space-time, which we will give in the next section. 

One of the powerful approaches to study a 2-dimensional fermion system is to 
utilize the bosonization technique,™ though it may be applicable only for charge- 
zero sectors of fermions, i.e for sectors of particle and anti-particle pairs. With this 
technique, one can expose hidden properties of such a fermion system by deriving the 



*' E-mail address: jsakamot@riko.shimane-u.ac.jp 



typeset using J>TpT|=X.sty <ver.0.8> 
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equivalent boson system. Since Coleman's discovery of the equivalence between the 
massive Thirring and the Sine-Gordon models, a* many people have tried to expand 
the freedom by introducing a non-Abelian group, by going into higher dimensional 
space-time, and so on. Formulations in most of these cases, however, turn out to be 
very complicated.™ 

In this article we apply the bosonization technique to the field theoretical GNF 
model with path integral quantization in order to investigate the physical structures 
of this model. As will be shown in the following, this is one of easier, but non-trivial, 
directions to expand the freedom of the system in the bosonization technique. We 
will study here the case of two species of fermions as the simplest case. Extension to 
the case of more than two species will be discussed elsewhere. We find in this model 
that for a specific value of the coupling constant one of the boson fields becomes 
free, while the other satisfies the Sine-Gordon equation. For that value of coupling 
constant, therefore, our two-species model is equivalent to the ordinary Sine-Gordon 
model of a single species. 

We use the metric convention in Minkowski space-time, g^ u = (— 1,+1) and 
e 01 = — eoi = 1. Gamma-matrices are given as 7 = —70 = ia x , 7 1 = 71 = a y , 
75 = 7o7i = 0"zj where a x , a y and a z are the Pauli matrices. 



§2. Model 



Our initial Lagrangian is given by 

£ = 4>a($ ~ m )V>a + lpb($ ~ m )^b + -jOJanOb, (2-1) 

where a and b denote the fermion species and 



In Q2-l| ), neither a- nor b- fermion interacts with itself directly. As mentioned in 
the previous section, one can derive the quantum mechanical GNF model from this 
Lagrangian. If we include a self-coupling term like (J2i=ab3i^) 2 instead of j a ^j^ in 
(2-1) above, which is one of the diagonal terms of SU (2), the model becomes simple, 
because we need only a single auxiliary vector field for such a case, unlike ( |2-3| ) below. 
We can suppose g > without loss of generality because we can change its sign by 
charge conjugation of one of fermion species, a or b. 

We introduce auxiliary vector fields and to rewrite the Lagrangian ( |2-1| ) 

as 

C = ^ a (@ - m)Va + M@ ~ m)A + \gA^ + - \gA^. (2-3) 

In 2-dimensional space-time, we can write vector fields and with scalar fields 
4> and x as 

Aft = e^ u d u (j) a + df+Xa, 

B^ = e^dTfo + d^Xb- (2-4) 
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Then we have 



C = + y^^T" + f fea - m)Va + (a - b) 

+ 9 -{d^ad^ b -d^ X ad^Xb)- 



(2-5) 



We transform the fermion variables as 



V>a -> = eX P <! ^f(-750a + Xa) [• V>a, 



-> V'f, = exp { —(-7506 + X&) }> V>b 



(2-6) 



to rewrite the Lagrangian as 



£' = - mex P{*5750a})V'a + ^&(# - wexp{i57 5 (?!> 6 })^6 



(2-7) 



Though transformation (2-6) causes the x to decouple from the other fields, it changes 
the path-integral measure of the fermion fields as 

J[ d^ a dip a = Y[ d^' a d^' a det | exp{-igj 5 c/) a } | , (2-8) 
and following Fujikawa's procedure & this determinant is calculated as 

det | exp{-igj 5 (f) a } |= exp ji J d 2 z|-(<90 a ) 2 j . (2-9) 

Then, the effective Lagrangian is given by 



£ c ff = i>JS ~ mexp{ig~f 5 (f> a })ip a + ij) b {$ - mexp{ig-f 5 (f> b })ip b 

2 

+ 9 -d^ a d^ b + |-{(^ a ) 2 + (dfo) 2 }, 

l 47T 

and the generating functional Z of the Green function is 

Z = J Y[ dipidvpidrfiiexpi J d 2 xC e R. 

i=a,b 



(2-10) 



(2-111 



Here and hereafter we write tjj for ip' . By the 'Wick' rotation we transform ourselves 
into Euclidean space-time from the Minkowski space-time. Then the generating 



functional ( 2-11 ) is rewritten as 

Ze = J IT dipidipid^i exp (— J d 2 



^(^-me is7A )4 + (a^b) 



+ 9 -d^ a d^ b + 9 —{ (d<t> a ) 2 + (dfo) 2 } 



(2-12) 
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§3. Perturbative expansion 

Before we expand Ze of (2-12) in power of m, we transform (0 o ,(fo) to (4>' a , <f>' b ) 

as 

<t> b =%4>a-<t>' b ) (34) 

to rewrite the Lagrangian (2-10) as 

C cS = ^ a {$ - me«*»)i + i) h (§ - me«^)^ 

+f(l + f)(W + f(l-f)(W- (3-2) 
Now, we make a perturbative expansion of in (2-12) in m: 

Ze = J JJ dipdipdcfr' exp ^— ^ d 2 x£ c ff (m = 0)^ 

x f^!^y d 2 x {^ ae ^{M b )^^ a + ^e^^-"^)^^) }" • (3-3) 



To obtain a boson model corresponding to our fermion model, we begin by calculating 
the first few terms of the expansion. It is seen that the n = 1 term vanishes due to 
the traceless property of the 7-matrices. The n = 2 term is given by 



Z E ] = J Yldi/fdipdtf/ exp(^- J d 2 x£ cS {m = 0) 

l J d 2 x ^e^^+^^^a + W^^-**)^^ 



+ ^Tr If^^^-^lT . ^igTOW^Dlj^ 
-2^Trr 1 e W5( ^ + ^ ) ^) (Tr^V^"^^} • (3-4) 

Here (• • •) denotes the vacuum expectation value of the time-ordered product, and 
Tr represents the trace with respect to space-time coordinates and 7-matrices. We 
have made use of the fact that the fermion propagator is given by 

It is easily seen that the third term of Eq. (3-4) vanishes due to the traceless property 
of 7-matrices, and that the first and second terms are equal. The first (or second) 
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term of Eq. (3-4) can be rewritten as 
8ir 2 J (x — y) 2 \ 

m? f o ,n 1 



8^J dxdy 



{x - y) s 



x exp |Wa Wa(2/)> " (^a 2 ) + " <tf >}■ (3-6) 

Here tr represents the 7-matrix trace. From Eq. (3-2) the boson propagators are 
given as 

i^'ai^Kiy)) = A 'ai x -y) = -A+ ln(x - y) V, 

(ti(x)4' b (y)) = A' b (x -y) = -A_ ln(x - y) V, (3-7) 

where \x is a small infrared cut-off mass, which will be set to zero after the calcula- 
tions, and we have put 

A± = — r. (3-8) 



Making use of the identity 



{x 



(3-9) 



we have 



z {2) = m 2 fl 2 {tf+gnW^+tf-g*)^)} 

E 2vr 2 



x J d 2 xd 2 yexpg 9 2 + ^{(^(xK(y))-(^ 2 )} 
xexp (j 5 2 - MixWM) - (<^ 2 )} 



2vr 2 



d 2 x(i 2 y ^cos s/Sg 2 /2 + g%{(f>' a (x) - <j)' a {y)} 

x ^cos s]Zg 2 /2-g-K{<t> h (x) - <&(y)}\ , (3-10) 

where we put0 

z = exp{(g 2 + gir)^ 2 ) + {g 2 - g^ 2 )}. (3-11) 
Setting m' 4 = 4m 2 fj, 2 z / 'it 2 , we can rewrite Eq. (3-10) as 

Zf = U[rn' 2 J d 2 xco S ^3g 2 /2 + g^' a cos ^3g 2 /2 - gir^J . (3-12) 



*' This term is, of course, divergent. We suppose that boson propagators are properly regularized 
to remove ultra-violet divergences. 
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This term is equal to the m /4 -order term of the perturbative expansion of the boson 
model given by 

+m' 2 cos I^SsVa + 9*4>' a \ cos {^3^/2 - gi<H\ - (3-13) 

One should note that the m /2 -order term in the expansion of the above model 
vanishes because of the super selection rule for the massless bosomfield in two- 
dimensional space-time, i.e. (exp^/3^) vanishes unless P = O.eP'EP In terms of 
the original boson fields, (j) a and 4>b, the Lagrangian above can be rewritten as 



+ m' 2 cos I ffi + f fa, + fc) I cos J ffi - f fa a - &) J . (3-14) 

In the Appendix we show that Ze of (2-12) is equal to that from (3-13) order 
by order in the perturbative expansion of m or m! . We conclude, therefore, that the 
fermion model of Eq. ( |2-1| ) is equivalent to the boson model of (3-13) or (3-14). 

§4. Discussion 

From the line of reasoning of the previous section, one might conclude that one 
can find the corresponding boson model by calculating only the first non-vanishing 
term of the perturbative expansion. But this is not the case for multi-species models 
like here. In fact, instead of (3-9), with the parameter a we can write 



{x - y) 



2 



^exp^A' a (x-y) + ^^A' b (x-y)Y (4-1) 



where a = 1/2 for (3-9). Then the Lagrangian of the corresponding boson model for 
this case seems to be given by 

£ = f(f+i)w 2 +!(f-i)w;) 2 
+ „ Aa j^j„j^j, „ 2) 

instead of (3-13). This ambiguity arises from the extra freedom due to the existence of 
two species of fermions, and does not exist in the original massive Thirring model. It 
is, however, easily seen that the Lagrangian above does not reproduce the expansion 
terms of (2-12) for higher orders than 0(m" 4 ) unless a = 1/2 (see the Appendix). 
The boson model (3T4), therefore, is a unique one which corresponds to the fermion 
model of fl2T|) . It would be interesting to ask to what order one should calculate to 
fix the corresponding boson model for the case of an n-species model. 



7 



Next we would like to refer to operator relations between fermion and boson 
fields. As usual, by comparing (3-14) with fl2T|) we may put 



V^M^fW) 2 , (4-3) 

Ipaln^a ^ud V (j) a , (4-4) 

and the same for b fields. These are the same as the original one-species case for the 
massive Thirring and Sine-Gordon models. We may also put 

m(tp a tp a + i>bA) m ' 2 cos ^3g 2 /2 + W>aj cos^3g 2 /2 - gir<j>' b ^ . (4-5) 

It is difficult to separate the right-hand side of the above expression into terms each of 
which contains only (f> a or <p b . We cannot, therefore, make correspondences between 
the fermion and boson fields for each species. It is hard to expect a correspondence 
like (4-5) from a naive extension of the original single-species case. 



We have stated nothing about value of g so far, except its sign. It seems that 
for 7T > g > 0, (f)' b becomes a ghost field, because the coefficient of its kinetic term in 
the Lagrangian (3-13) is negative for such values. From the viewpoint of the path 
integral formulation, we see that the quantity J d(j)' b exp{—S(m = 0)} is ill-defined 
even if we are in Euclidean space-time. To investigate a two-body composite system, 
use of the Bethe-Salpeter (BS) equation is often made. It is known that there is an 
abnormal solution for the BS equation for most cases. We conjecture here that the 
^-field would correspond to such an abnormal solution for the particle-anti-particle 
composite state. 

We find from (3T3) that <p' b becomes a free field for g = 2ir/3, though it is in 
the ghost region mentioned above, while <j)' a becomes a Sine-Gordon field. Therefore, 
for such a specific value of g, our system is equivalent to an ordinary Sine-Gordon 
model with a single species. A similar fact is seen for some SU{2) extension of the 
model. As mentioned above, if we take (J2i= a b 3i) 2 > which is one of the diagonal 
terms of SU(2), instead of j a ^jb i n fl2T| ), we only need a single auxiliary bose field, 
unlike (|2-3| ). The value of g, however, is not concerned with this case. 

If we rescale the boson fields as 



ff*/2n + g/2<t/ a , (4-6) 



- <t>'l = - g/201 (4-7) 

then we have from (3-13) 



r , 1 /^ 5 J 1 //^2 , 2m ' I / 37Tg + 2tt 2 

^=2^ d ^) +2^) + cos W g + ~ ^ ( cos 



For g = we have 

,/2 




C = -(dc2> + ) 2 + -(<9<2>_) 2 + — (cos^ : <Z> + + cosv / 4^>-}, (4-9) 

2 2 7T ^ J 
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where we have set 

*± = ^(#±tf). (4-10) 
This is the well-known boson model which corresponds to the free fermion model. 
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Appendix A 

The n-th Order Terms of the Perturbative Expansion 



We now prove that the generating functional Ze of (2-12) is equal to that of 
(3-13) order by order in the expansion with respect to m or m! . The n-th order term 
of Ze is given as 

4" } = / II dxd V ({^ a e l9 ^a + ? 6 e**^ 6 } n ) . (A-l) 

As is mentioned in the text, the odd order terms of the above expansion vanish. 
Then we set n = 2k to obtain 

7 (2k) _ m 2k ^ (2fc)! 



(2k)\^(2r)\(2k-2r)\ 
x ( / dx{a -> b)Y ) . (A-2) 



Due to the super selection rule for massless boson fields, the number of +<p a (+<&>) 
must be equal to that of — 4> a (—4>b) in exponentials of the above expression, and 
thus we have 



■r- 



-(H)H(fc-r)!} 2 



k— r\ 



m 2k 



r=0 v LV rs i=l 

I r 

X 



4 = 1 " 
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f[ ^b( x j)—^M x j)^b(yj)^r^Myj) 



\j=r+l 



x ( exp ig 



x 2k k 



XI {^a(xj) - 4> a {yi)} + X ~ Mvj)} 

i=l j=r+l 



1 



2W ^( r !)2 { ( fc _ r)l} s 



1=1 

n[>j=l( x t ~ x j) 2 (Vi — Vj) 2 Ui>j=r+l( x i ~ x j) 2 (Vi ~ Vj) 2 

Ulj=i( x i - yj) 2 

x ( exp i<7 



IIij'=r+l( a 'i 



X {0a(a*) - MVi)} + E - MVj)} 

i=l j=r+l 

On the other hand, from the boson model (3-13) we have 

( n ) m' 2n f n I n \ 

Z E(B) = J Tl dx i ( J[cOS (3 + <p' a (xi) COS f3-(p' b {xi) j 

-i / i s 2n - n In \ 

=Mt) /n^(nf M,i) +^ M( l ( 

where we have put 



(A-3) 



(A-4) 

P± = V 3 ^ 2 /2 ± 5vr. (A-5) 
We can set n = 2k, because only terms in the above expression containing the same 
numbers of +4>' a and —<j>' a in exponential do not vanish. We have 



/ 2k 

n{< 



,i/3+(p' a (xi) + e -i/3 +( f> 



k\k\ 



expi^+^K^O-to)}), (A-6) 



where we have rearranged arguments {xi, X2, • • • , X2fc} — ► {xi, X2, ■ ■ ■ , xjt, yi, • • • , y^} 
so that the Xj and are arguments of the an d { — 'P'a}-, respectively. By this 

rearrangement, the last factor of (A-4) becomes 



E 



k\ 



r = \r\(k — r)\) \ 
Then we have 

k 



( exp i(3- 



- ^ 6 (w)} - E Mte) - ^(W)} 

i=l i=r+l 



4?2) = (™') 4fc E (^3751) / R dxdy ( expi/3 + E {te) - &(w)} 



x ( exp 



EWW-to)}- E - to)} 

i=l i=r+l 



(A-7) 
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Making use of the identities 

P 2 + A' a (x -y) + 2 _A' b (x -y) = g 2 '(M^Mv)) ~ H* - y)V , 

and 



f3 2 + A' a {x -y)- (3iA' b (x -y) = g' {4> a (x)<j> b {y)) , 



(A-8) 
(A-9) 



we obtain 



7 {2k) _ 



m\ 2k " f 1 



r=0 



n[>j=i(^ - gj) (yt - rii>j=r+i(^ - xj) (vi - vj) 



x exp g 



^2 {A a {xi - Xi) + A a (yi - yj)} 
i>j=i 



+ {A a (xi- xj) + A a (y i -y j )} 

i>j=r+l 

r k 

- ^2 A a (xi - yj) - ^a{xi-yj) 

i,j=l i,j=r+l 
r k 

+ i A ab{xi - Xj) + A ab (yi - yj) - A ab (xi - yj) - A ab (xj - yi)} 

i=l j=r+l 

(A- 10) 

where we put 

A a ( x - y )= (Mx)Mv)) - (Mx) 2 ) = (Mx)Mv)} - (M*) 2 ), (a-ii) 



Aab(x ~ y)= ((p a (x)(t>b(y)) ~ ((/>a(x)(/) b (x)). 



(A-12) 



Then we find Eq. (A- 10) is equal to (A-3), i.e. Z^ k ' = Zgfyy 

In the above argument it is essential that Eqs. ( |AS| ) and ( |A-9| ) hold. If we take 
(4-2) as the corresponding boson model, then putting 



2a( 5 2 + gir) + 5 2 /2, 



(A-13) 
(A- 14) 



i2(l-a)(g*-g7T)+gy2 
instead of (3 + and /3_ respectively, we have 

n\A' a {x -y) + V-A b (x -y) = g 2 (M^)Mv)) ~ ~ vf V \ (A-15) 
V 2 + A' a (x -y)- rlA' b (x - y) = g 2 (Mx)My)) ~ (2" - 1) m(x - y) V, (A- 16) 
instead of (|A-8|) and ( |A-9|) . Due to the last term of ( |A-16|) , Z^ k ' is not equal to 

^E(B) UIHess a = 1/2- 
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